Construction of copy number variation landscape and characterization of associated genes in a Bangladeshi cohort of neurodevelopmental disorders.

Introduction: Copy number variations (CNVs) play a critical role in the pathogenesis of neurodevelopmental disorders (NDD) among children. In this study, we aim to identify clinically relevant CNVs, genes and their phenotypic characteristics in an ethnically underrepresented homogenous population of Bangladesh. Methods: We have conducted chromosomal microarray analysis (CMA) for 212 NDD patients with male to female ratio of 2.2:1.0 to identify rare CNVs. To identify candidate genes within the rare CNVs, gene constraint metrics [i.e., "Critical-Exon Genes (CEGs)"] were applied to the population data. Autism Diagnostic Observation Schedule-Second Edition (ADOS-2) was followed in a subset of 95 NDD patients to assess the severity of autism and all statistical tests were performed using the R package. Results: Of all the samples assayed, 12.26% (26/212) and 57.08% (121/212) patients carried pathogenic and variant of uncertain significance (VOUS) CNVs, respectively. While 2.83% (6/212) patients' pathogenic CNVs were found to be located in the subtelomeric regions. Further burden test identified females are significant carriers of pathogenic CNVs compared to males (OR = 4.2; p = 0.0007). We have observed an increased number of Loss of heterozygosity (LOH) within cases with 23.85% (26/109) consanguineous parents. Our analyses on imprinting genes show, 36 LOH variants disrupting 69 unique imprinted genes and classified these variants as VOUS. ADOS-2 subset shows severe social communication deficit (p = 0.014) and overall ASD symptoms severity (p = 0.026) among the patients carrying duplication CNV compared to the CNV negative group. Candidate gene analysis identified 153 unique CEGs in pathogenic CNVs and 31 in VOUS. Of the unique genes, 18 genes were found to be in smaller (<1 MB) focal CNVs in our NDD cohort and we identified PSMC3 gene as a strong candidate gene for Autism Spectrum Disorder (ASD). Moreover, we hypothesized that KMT2B gene duplication might be associated with intellectual disability. Conclusion: Our results show the utility of CMA for precise genetic diagnosis and its integration into the diagnosis, therapy and management of NDD patients.

Analyzing single cell transcriptome data from severe COVID-19 patients.

We describe the protocol for identifying COVID-19 severity specific cell types and their regulatory marker genes using single-cell transcriptomics data. We construct COVID-19 comorbid disease-associated gene list using multiple databases and literature resources. Next, we identify specific cell type where comorbid genes are upregulated. We further characterize the identified cell type using gene enrichment analysis. We detect upregulation of marker gene restricted to severe COVID-19 cell type and validate our findings using in silico, in vivo, and in vitro cellular models. For complete details on the use and execution of this protocol, please refer to Nassir et al. (2021b).

Mutational Landscape of Autism Spectrum Disorder Brain Tissue.

Rare post-zygotic mutations in the brain are now known to contribute to several neurodevelopmental disorders, including autism spectrum disorder (ASD). However, due to the limited availability of brain tissue, most studies rely on estimates of mosaicism from peripheral samples. In this study, we undertook whole exome sequencing on brain tissue from 26 ASD brain donors from the Harvard Brain Tissue Resource Center (HBTRC) and ascertained the presence of post-zygotic and germline mutations categorized as pathological, including those impacting known ASD-implicated genes. Although quantification did not reveal enrichment for post-zygotic mutations compared with the controls (n = 15), a small number of pathogenic, potentially ASD-implicated mutations were identified, notably in TRAK1 and CLSTN3. Furthermore, germline mutations were identified in the same tissue samples in several key ASD genes, including PTEN, SC1A, CDH13, and CACNA1C. The establishment of tissue resources that are available to the scientific community will facilitate the discovery of new mutations for ASD and other neurodevelopmental disorders.

Single-cell transcriptome identifies FCGR3B upregulated subtype of alveolar macrophages in patients with critical COVID-19.

Understanding host cell heterogeneity is critical for unraveling disease mechanism. Utilizing large-scale single-cell transcriptomics, we analyzed multiple tissue specimens from patients with life-threatening COVID-19 pneumonia, compared with healthy controls. We identified a subtype of monocyte-derived alveolar macrophages (MoAMs) where genes associated with severe COVID-19 comorbidities are significantly upregulated in bronchoalveolar lavage fluid of critical cases. FCGR3B consistently demarcated MoAM subset in different samples from severe COVID-19 cohorts and in CCL3L1-upregulated cells from nasopharyngeal swabs. In silico findings were validated by upregulation of FCGR3B in nasopharyngeal swabs of severe ICU COVID-19 cases, particularly in older patients and those with comorbidities. Additional lines of evidence from transcriptomic data and in vivo of severe COVID-19 cases suggest that FCGR3B may identify a specific subtype of MoAM in patients with severe COVID-19 that may present a novel biomarker for screening and prognosis, as well as a potential therapeutic target.

Single-cell transcriptomics trajectory and molecular convergence of clinically relevant mutations in Brugada syndrome.

Brugada syndrome (BrS) is a rare, inherited arrhythmia with high risk of sudden cardiac death. To evaluate the molecular convergence of clinically relevant mutations and to identify developmental cardiac cell types that are associated with BrS etiology, we collected 733 mutations represented by 16 sodium, calcium, potassium channels, and regulatory and structural genes related to BrS. Among the clinically relevant mutations, 266 are unique singletons and 88 mutations are recurrent. We observed an over-representation of clinically relevant mutations (∼80%) in SCN5A gene and also identified several candidate genes, including GPD1L, TRPM4, and SCN10A. Furthermore, protein domain enrichment analysis revealed that a large proportion of the mutations impacted ion transport domains in multiple genes, including SCN5A, TRPM4, and SCN10A. A comparative protein domain analysis of SCN5A further established a significant (P = 0.04) enrichment of clinically relevant mutations within ion transport domain, including a significant (P = 0.02) mutation hotspot within 1321-1380 residue. The enrichment of clinically relevant mutations within SCN5A ion transport domain is stronger (P = 0.00003) among early onset of BrS. Our spatiotemporal cellular heart developmental (prenatal to adult) trajectory analysis applying single-cell transcriptome identified the most frequently BrS-mutated genes (SCN5A and GPD1L) are significantly upregulated in the prenatal cardiomyocytes. A more restrictive cellular expression trajectory is prominent in the adult heart ventricular cardiomyocytes compared to prenatal. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.NEW & NOTEWORTHY Brugada syndrome is a rare inherited arrhythmia with high risk of sudden cardiac death. We present the findings for a molecular convergence of clinically relevant mutations and identification of a single-cell transcriptome-derived cardiac cell types that are associated with the etiology of BrS. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.

Long-Read Sequencing Improves the Detection of Structural Variations Impacting Complex Non-Coding Elements of the Genome.

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value < 1.1020 × 10-57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1-100 kb size range (p-value < 9.080 × 10-67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value < 9.02 × 10-9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome.

Whole exome sequencing uncovered highly penetrant recessive mutations for a spectrum of rare genetic pediatric diseases in Bangladesh.

Collectively, rare genetic diseases affect a significant number of individuals worldwide. In this study, we have conducted whole-exome sequencing (WES) and identified underlying pathogenic or likely pathogenic variants in five children with rare genetic diseases. We present evidence for disease-causing autosomal recessive variants in a range of disease-associated genes such as DHH-associated 46,XY gonadal dysgenesis (GD) or 46,XY sex reversal 7, GNPTAB-associated mucolipidosis II alpha/beta (ML II), BBS1-associated Bardet-Biedl Syndrome (BBS), SURF1-associated Leigh Syndrome (LS) and AP4B1-associated spastic paraplegia-47 (SPG47) in unrelated affected members from Bangladesh. Our analysis pipeline detected three homozygous mutations, including a novel c. 863 G > C (p.Pro288Arg) variant in DHH, and two compound heterozygous variants, including two novel variants: c.2972dupT (p.Met991Ilefs*) in GNPTAB and c.229 G > C (p.Gly77Arg) in SURF1. All mutations were validated by Sanger sequencing. Collectively, this study adds to the genetic heterogeneity of rare genetic diseases and is the first report elucidating the genetic profile of (consanguineous and nonconsanguineous) rare genetic diseases in the Bangladesh population.